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For the quantum Ising model with ferromagnetic random couphngs Jij > and random trans- 
verse fields hi > ai zero temperature in finite dimensions d > 1, we consider the lowest-order 
contributions in perturbation theory in {Ji,j/hi) to obtain some information on the statistics of 
various observables in the disordered phase. We find that the two-point correlation scales as : 

In C(r) ~ ^ r^u, where S^typ is the typical correlation length, it is a random variable, and cj 

^ ■ coincides with the droplet exponent ujdp{D = d — 1) of the Directed Polymer with D = [d — 1) 

_ transverse directions. Our main conclusions are (i) whenever cj > 0, the quantum model is governed 

by an Infinite-Disorder fixed point : there are two distinct correlation length exponents related 

■ by Vtyp = (1 — uj)vav the distribution of the local susceptibility Xioc presents the power-law tail 
\ P{Xioc) ~ ^/Xioc where vanishes as ^au", so that the averaged local susceptibility diverges in 

a finite neighborhood < /i < 1 before criticality (Griffiths phase) ; the dynamical exponent z 

■ diverges near criticality as z = dj [i ~ (ii) in dimensions d < 3, any infinitesimal disorder 
flows towards this Infinite-Disorder fixed point with > (for instance a;(d = 2) = 1/3 and 
Ltj(d = 3) ~ 0.24) (iii) in finite dimensions d > 3, a finite disorder strength is necessary to flow 
towards the Infinite-Disorder fixed point with a;(d) > (for instance a;(d = 4) ~ 0.19), whereas a 
Finite-Disorder fixed point remains possible for a small enough disorder strength. For the Cayley 
tree of effective dimension d = oq where cj = 0, we discuss the similarities and differences with the 

C/3 ■ case of finite dimensions. 
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INTRODUCTION 



The quantum Ising model 

o 

involving a ferromagnetic coupling J > between nearest neighbors < i,j >, and a transverse field h on each site, 
■ is the basic model to study quantum phase transitions at zero-temperature This model on an hypercubic lattice 
<*' . in dimension d is well understood via the equivalence with the classical Ising model in dimension = d + I, 

i.e. the time plays the role of an extra space-dimension [T], and the dynamical exponent is Zp^re = 1- The lower- 
critical dimension df""^" = 1 and the upper critical dimension c?^/""* — 4 for the classical model yields di = and 



• du — S for the quantum model. Above the upper critical dimension d > du = 3, the mean-field critical exponents 
{aMF — 0',Pmf — 1/2; 7mf — ^',i^mf = 5) are exact. In dimension rf = 1, the critical exponents are exactly known 
from the d'^'"*'' — 2 classical Ising model. 

Let us now consider the disordered model, where the nearest-neighbor couplings Jij > and the transverse-fields 
hi > are independent random variables drawn with two distributions 'k coupling {J) and TTfieidih) 



^ = - E ^M<'^i-E^^< (2) 

<ij> i 

In dimension d = 1, exact results for a large number of observables have been obtained by Daniel Fisher [2| via the 
asymptotically exact strong disorder renormalization procedure (for a review, see 0]). In particular, the transition 
is governed by an Infinite-Disorder fixed point and presents unconventional scaling laws with respect to the pure 
case. In dimension d > 1, the strong disorder renormalization procedure can still be defined. It cannot be solved 
analytically, because the topology of the lattice changes upon renormalization, but it has been studied numerically 
with the conclusion that the transition is also governed by an Infinite-Disorder fixed point in dimensions c? = 2, 3, 4 [j- 
[T^ . These numerical renormalization results are in agreement with the results of independent quantum Monte-Carlo 
in d = 2 15, 16]. 

In dimension d = 1, one has the important additional property that any infinitesimal disorder fiows towards the 
Infinite-Disorder fixed point [5, 's^. To try to answer this question in dimension d > 1, the first natural step is to 
determine the relevance of disorder at the transition of the pure model via the use of the Harris criterion [iTj or the 
inequality 1/ > 2/d of Chayes et al |18|. Below the upper-critical dimension d < d„ = 3 of the pure quantum model, 
there exists a single correlation exponent Vpure, and the conclusion is that disorder is relevant for Vpure < 2/d, and 
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irrelevant for Vpure > 2/(i. Here, the conclusion is thus that disorder is relevant in dimension d = 1 where Vpureid = 
1) = v^^^^'^{(r^'"''' = 2) = 1< 2/(i = 2, in dimension d = 2 where Vpure{d = 2) = v^^^^'^ {d"^"'''' = 3) ~ 0.63.. < 2/d = 1, 
and in dimension d — ?> where Vpure{d = 3) = v^^^^ {d'^'-"-^'^ — 4) = 1/2 < 2/d — 2/3. Above the upper-critical 
dimension d > = 3 of the pure model, there exists two correlation length exponents : the usual one v^^^ = 1/2 
defined from the two-point correlation, and the exponent vi that governs finite-size scaling, where vf"-^^ = 2/d'^'°*'' 
for the classical Ising model in dimension d'^'"'^" > d^'""* = 4 [l^, and where j^?"''"*"™ — 3/{2d) for the quantum 
Ising model in dimension d > du = 3 2Q-] (because the quantum model in dimension d is equivalent to a (d -f- 1) 
classical model in a special cylindrical geometry, see pd | and references therein). So the use of the Harris criterion is 
more involved, since the two values f and fi can give different answers [2ll |. In our present case, the Harris criterion 
based on v = 1/2 yields that disorder should be irrelevant for d > 4 [H, whereas the Harris criterion based on 

i/i = 3/ (2d) leads to the conclusion that disorder should be always relevant ! In addition, these arguments based on the 
weak-disorder Harris criterion have been questioned since rare regions can play an essential role in these quantum 
disordered models, where the disorder is actually 'infinitely' correlated along the time-direction. The effects of rare 
regions are discussed in detail in the review [22], with the conclusion that the important parameter is the effective 
dimensionality d^jj of rare regions (in our present quantum model where the disorder is actually 'infinitely' correlated 
along the time-direction, the dimensionality of rare regions is djfjf = 1 p2j). This dimensionality djiji of rare regions 
should be compared to the lower critical dimension df'^'^'' sufficient to obtain ordering (in our present quantum model 
where the order is a ferromagnetic magnetization, the lower critical dimension is df""^^ — The random quantum 
Ising model thus corresponds to the case dRR = df""^^, i.e. to the so-called class B in the classification described in 
[23 |. where rare regions dominate the critical behavior and induce an unusual activated scaling. (The conventional 
power- law scaling is expected to hold for the so-called class A corresponding to dim < df"-^^, where rare regions 
cannot undergo the phase transition by themselves, whereas the so-called class C corresponds to the case dRR > df"'''^ 
where rare regions can order by themselves at different values of the order parameter) . In conclusion, this argument 
based on rare regions suggests that the unusual activated scaling found by Daniel Fisher ^ in dimension d = 1 via 
the asymptotically exact strong disorder renormalization procedure should exist in any finite dimension d and even 
on the Cayley tree (d — oo). 

To elucidate what happens in high finite dimension d, it is natural to consider the limit of infinite dimension 
d = oo. Two opposite conclusions have appeared recently : on one hand, an Infinite-Disorder fixed point has been 
found numerically for the Erdos-Renyi random graph in [l3], but on the other hand, the solution on the Cayley tree 
obtained via some approximations within the quantum cavity framework [23l - [25| has been interpreted as a 'Finite- 
Disorder critical point' (with however strong inhomogeneities). 

In this paper, our aim is to clarify the possibility of an 'Infinite Disorder' fixed point as a function of the dimension 
d via simple scaling arguments. We focus on the disordered phase of the random quantum Ising model Eq. [51 and we 
consider the leading contributions in perturbation theory in the disordered variables Jij/hi for various observables. 
We conclude that the properties of the quantum model in its disordered phase are determined by the problem of the 
Directed Polymer in a random medium with D = {d — 1) transverse directions. We analyze the various consequences 
of this correspondence. Since the solution on the Cayley tree obtained via the approximated quantum cavity method 
(23l-|25j is also related to the problem of the Directed Polymer on a tree, we rediscuss this case to emphasize the 
similarities and differences between finite dimensions and infinite dimension. 

The paper is organized as follows. In section|lTl we compute the magnetization of a center-site when an infinitesimal 
magnetization rriext — >■ is imposed on the surface of a finite sample, via a perturbative expansion deep in the 
disordered phase, and we explain the link with the scaling properties of the Directed Polymer model. In section Hill 
we analyze the consequences for the two-point correlation, and obtain that there are two distinct correlation length 
exponents Utyp and Vav for the typical and averaged correlation lengths. In section HVl we discuss the power-law tail 
of the distribution P(xioc) ^ ^/Xioc^ of the local susceptibility x;oc, where /i vanishes near criticality. In section IVl 
we analyze the statistics of the gap and find that the dynamical exponent diverges near criticality as z = d/fj,. Finally 
in section I VII we consider the case of the Cayley tree to compare with the case of finite dimensions d. Section IVIII 
summarizes our conclusions. Appendix |^ contains a reminder on properties of the Directed Polymer model that are 
used in the text. In Appendix |B1 we discuss the analogy with Anderson localization, where the droplet exponent of 
the Directed Polymer also appears in the localized phase. 
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II. MAGNETIZATION OF A CENTER SITE 
A. Perturbative expansion deep in the disordered phase 

When the ferromagnetic couphngs vanish Jij = 0, all spins are completely independent, there is no order even at 
the distance of one lattice spacing. Deep in the disordered phase of the quantum model, where the correlation length 
remains short, it is thus natural to consider this uncorrelated state as the reference state, and to study a perturbative 
expansion in the couplings Jij. 

In this section, we consider the magnetization toq =< cr^(O) > of a center site in the ground state of a finite 
sample of volume (2L)'^ where an infinitesimal magnetization uiext — >■ is imposed on the surface S of the sample. At 
first order in the couplings J^^-, the magnetization ttiq is related to the magnetizations uij of its K = {2d) neighbors 
on the hypercubic lattice 

2d J 



= 1 



If we iterate up to the boundaries to reach the imposed infinitesimal magnetization niext — ^ 0, we obtain a sum over 
the paths P going from the center site to the surface S 

PathP:6^S 

where the weight W{P) of each path P is the product of the factors ) along the path 




W{P)^ ^ (5) 



B. Equivalence with the Directed Polymer model 

In the whole disordered phase, the magnetization mo of the center site should typically decay exponentially with 
respect to the distance L to the boundary 

In — + ... (6) 

where £,typ represents the typical correlation length of the disordered phase. 

Deep in the disordered phase where we have written the perturbative expansion of Eq. 21 we thus expect that the 
weight W{P) of each path P will typically decay exponentially with the length 1{P) of the path. As a consequence, 
the sum over paths in Eq. |4] will be dominated by directed paths going from the center site to the surface at distance 
L. So the scaling properties will be the same as for the Directed Polymer in a random medium with D = (d — 1) 
transverse directions (see Appendix \X\ for a reminder on the Directed Polymer model). In particular, whenever the 
Directed Polymer is in its localized phase characterized by the droplet exponent u = ujdp{D) (see section lA II in 
Appendix), the center-site magnetization will presents the scaling 

In^ ^ ~^ + L-u (7) 

where m is a random variable of order 0(1). 

Deep in the disordered phase, the mapping to the Directed Polymer model should be quantitative, i.e. toq of Eq. 
|4] corresponds to the partition function of Eq. I All at temperature T = I with effective random energies £<,(«) on sites 
i determined by the transverse fields hi 

e,{i)= In h{i) (8) 
and with effective random energies eb{i,j) on bonds determined by the ferromagnetic couplings J^j 

e&(«7j) = - In Jij (9) 
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From the dimensionality D = d — I and the disorder distribution of this effective Directed Polymer model, one can 
determine whether the Directed Polymer is in its localized phase or not. From the reminder of section lA 3[ we 
conclude : 

(i) In dimension d < 3 for the quantum model, there is no delocalized phase for the associated Directed Polymer 
model, so that any initial disorder should correspond to the localized phase of the Directed Polymer characterized by 
the droplet exponent uj{d) = ujdp{D = d — 1). In particular, for the quantum model in dimension d = 2, we conclude 
that the corresponding exponent should take the explicit value (Eq. IA3|) 

uj{d = 2)=u;DpiD = l) = ^ (10) 

whereas in d = 3, it should be of order uj{d = 3) = ujdp{D = 2) ~ 0.24 (Eq. IA5[) . 

(ii) In dimension d > 3 for the quantum model, the associated Directed Polymer model can be in two different 
phases, depending on the strength of the disorder. For instance, if the effective energies of Eq. [S] and [S] are Gaussian 
(this corresponds to a log-normal distribution for the random fields and for the random bonds of the quantum models), 
the Directed Polymer will be in its localized phase only if the variances of the Gaussians are above some threshold. 
Then the quantum model will flow towards an Infinite-Disorder fixed point characterized by a;(d) = ijJdp{D — d — \) 
(see Eq. IA5I for numerical estimates in dimension d < 8). Below some disorder threshold, the associated Directed 
Polymer will be in its delocalized phase, and a Finite-Disorder fixed point for the quantum model is possible. Note 
however that in the quantum model, the most frequently used disorder distribution is a box distribution for the 
random fields 

7rso.(M = ^^(0</i<A) (11) 

Then the annealed partition function for the effective Directed Polymer model does not exist ( the inverse moment 
diverges J dh'iTBox{h)/h = -\-oo as a consequence of the finite weight at zero-field TTsoxih — 0) > ). Then there is no 
delocalized phase (see section 1X2)) and the Directed Polymer can only be in its localized phase. 

C. Extension to the whole disordered phase 

Up to now, we have discussed what happens deep in the disordered phase where the typical correlation length S^typ 
introduced in Eq. [7] is very small. Let us now discuss what can be extended to the whole disordered phase and 
what cannot. As the typical correlation length ^typ grows, the 'directed constraint' on paths will become less and less 
effective, higher orders of perturbation should be taken into account, and so on. As a consequence, no quantitative 
conclusion can be obtained for the typical correlation length ^typ, and in particular, its divergence at the critical point 
of the quantum model (where S represents some appropriate distance to the quantum transition) 

^typ cx S-'^^y- (12) 

is not determined by the physics of the Directed Polymer. 

However, we expect that the fluctuation exponent uj obtained previously deep in the disordered phase should remain 
the same in the whole disordered phase of the quantum model. As a comparison, we describe in Appendix [B] the case 
of Anderson localization, where the droplet exponent of the Directed Polymer also appears in the localized phase for 
completely similar reasons, whereas the position of Anderson transition and the divergence of the localization length 
cannot be predicted from the Directed Polymer model. 

Let us now discuss the special case of the dimension d = 1, where the validity of the perturbative expansion can be 
analyzed precisely via a direct comparison with exact result. 

D. Case d=l 

As recalled in the introduction, the quantum model in dimension d = 1 has been exactly solved by D. Fisher Q 
via the strong disorder renormalization procedure. Here to discuss the domain of validity of the above perturbative 
expansion, it is however simpler to consider not the magnetization of a center site, but instead the surface magnetiza- 
tion rriQ^^^ at the boundary 0, when an external magnetization mj^ — rriext is imposed at the other boundary of the 
finite chain of length L. 
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1. Surface magnetization at lowest order in perturbation theory in the couplings 



Within the perturbation theory described above, there exists a single directed path going from site to L on the 
chain, so that the surface magnetization mi^ reads at lowest order in perturbation theory 



sur f 
^0 



rUext 



L-1 

n 

i=0 



rriext exp 



L-l 



^(InJ, 



In hi) 



(13) 



The logarithm of the magnetization is given by a sum of random variables and thus displays the Central Limit scaling 

(14) 



/ sur f ^ 

In 



L(ln Ji^i+i —Inhi) + L^u 



i.e. this perturbative analysis yields that the fluctuation exponent is simply the random walk exponent (that corre- 
sponds to the degenerate case of a Directed Polymer with D = transverse dimensions) 



uj{d = 1) = UJDP{D = 0) = ^ 



(15) 



Comparison with the exact expression for the surface magnetization 



Let us now compare with the exact expression of the surface magnetization when one imposes = 1 [2 61] , that 
can be obtained from a free-fermion representation 



sur f 



L-1 i 

En 

i=0 j=0 



J. 



^1/2 



(16) 



At lowest order in the couplings Jjj+i, the sum is dominated by the term i = L — 1 containing the product of L terms 
and one recovers the perturbative expression of Eq. 1131 The statistical properties of the exact expression of Eq. [16] 
in both phases and at criticality have been discussed in detail in [26l - [28j . The important point is that the fluctuation 
exponent of Eq. [15] present in the perturbation theory remains the exact exponent in the whole disordered phase, as 
well as the Gaussian distribution of the variable u in Eq. 1141 Only the numerical prefactors of the linear term and of 
the fluctuation terms are renormalized with respect to the perturbative results. We note t hat here the pert urbative 
expression of Eq. [14] actually predicts also correctly the exact critical point determined by (In J^^i+i — \nhi) = and 
the typical correlation length exponent i^typ = 1, but we believe that these two last properties are peculiar to the 
dimension d — 1, and that in higher dimension d > I, the position of the critical point and the exponent vtyp cannot 
be obtained from the Directed Polymer problem (see the discussion in section III C[) 



E. Case of the Cayley tree 

Let us now consider the case where the quantum model is defined on a tree of coordination number {K + 1). The 
computation of the center site magnetization toq in terms of the infinitesimal external magnetization ruext imposed 
at generation L corresponds at leading order in perturbation theory to the Directed Polymer model defined on a tree, 
with the disorder parameters given by Eqs [8] and [9] : this corresponds exactly to the 'linearized equations' at zero 
temperature obtained after some approximations within the quantum cavity method ( Eq 5 of Ref [23| , Eq 9 of Ref 
[2^ or Eq 24 of [1^). We refer to Ref [l^ for the comparison of various levels of approximations. From the exact 
solution of the Directed Polymer on a tree [1^ , the droplet exponent vanishes (Eq IA4I in Appendix) 

= (17) 

whereas in all low finite dimensions D = d — 1 < 7 that have been tested numerically, the droplet exponent remains 
positive u!{d) = ujdp{D) > (Eq. IA5p . Since the two cases a; > and a; = induce very different scaling properties, 
we have chosen for clarity to focus on the case of finite d with uj{d) > in the following sections, and to discuss 
separately the case of the Cayley tree in section [Vl] 
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III. STATISTICS OF THE TWO-POINT CORRELATION WHEN uj{d) > 
A. Perturbative expansion deep in the disordered phase of the quantum model 

Deep in the disordered phase of the quantum model, the two-point correlation function 

C(f) ^< a£a£ > (18) 
has for leading contribution in perturbation theory (as in Eq. |4|) 

PathP:6-fr 

= ( n I (19) 

where the sum is over the paths P going from the point to the point r. In the disordered phase of the quantum 
model where the correlation decays exponentially, we expect again that directed paths dominate the sum, so that the 
scaling of the Directed Polymer model should again appear as Eq. [7] 

In C(f) ~ --!— + r'^u (20) 

r->+oo i^f^yp 

This generalizes the known solution in o? = 1 Q where uj{d = 1) ~ uj£)p{D = 0) = 1/2. For the quantum model in 
dimension d = 2, we thus expect uj{d = 2) = ujdp{D = 1) = 1/3 (Eq. [TU| . 

Besides the droplet exponent uj already discussed, it is now important to characterize the tail of the distribution 
P{u) of the random variable u 

InP(u) ~ -m" (21) 

From the results recalled in the Appendix, one has the simple relation (see the discussion around Eq. lASp 

v{d)=Tjnp{D = d-l) = - ,n , ^, = . \. (22) 

1 — UJDP{1J ^ d~ 1) 1 — uj{d) 

that generalizes the exponent r]{d = 1) = rj£ip{D = 0) = 2 of the Gaussian distribution in d = 1 (Eq. [H)) . 



B. Averaged correlation 

From the distribution of the logarithm of the two-point correlation of Eq. 1201 one may estimate the contribution 
of rare events to the averaged correlation : the correlation C(r) can be of order 0{1) if the random variable u takes 
the anomalously large value 

ui ~ +- (23) 

Z,typ 

and this event occurs with the small probability P{ui) given by 

j,(l-w)r) 

InP(wi) cx -u? ~ -75 (24) 

Since (1 — uj)ri = 1 fEq [22)1 . this probability of rare events where C(r) ~ 1 decays exponentially with the distance 
r, with a characteristic length S^^y^ that diverges more rapidly than the typical correlation length S^typ whenever the 
tail exponent satisfies ry > 1, i.e. whenever the droplet exponent does not vanish w > (Eq lASp . This means that the 
averaged correlation function Cai,(r) — C{r) will be dominated by these rare events and will decay exponentially 

lnCa„(r)~lnP(ui)~--^ (25) 

sow 
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where the averaged correlation length (,av diverges as 

1 

£,aV OC £,fyp — ^fyp (26) 

i.e. that there are two distinct correlation length exponents related by 

l^av = l^typV = (27) 
1 — OJ 

In summary, we expect that in finite dimension d, where the droplet exponent ujdp of the Directed Polymer does 
not vanish (see Eq. IA5I) . the disordered phase of the quantum model is characterized by two distinct correlation 
length exponents related by Eq. I27[ that generalize the well-known results vtyp — 1 and i^av = 2 in d = 1 . 

C. Finite-size scaling of the typical correlation 

When there exists a single correlation length, this length governs all finite-size properties near criticality. When the 
typical and averaged correlation lengths diverge differently, one expects that it is the averaged correlation length that 
governs finite-size effects. It is thus natural to write the following finite-size scaling form for the typical correlation Q 

hTcW ^ -^^"^^(7^) (28) 

ujc is the exponent governing the decay exactly at criticality where the scaling function reduces to some positive 
constant F{0) > 0. The matching with the scaling of Eq. ^U\m the disordered phase imposes the asymptotic behavior 



F{x) (X x^''^" (29) 



leading to 



InC(r) ~ ^ (30) 

i.e. by consistency with Eq. [20l the typical correlation length should be ^typ ~ ^Iv'^"- The comparison with Eq. [26] 
leads to the conclusion 

Wc = w (31) 
that generalizes the known result uJc = uj = l/2md=l In dimension d = 2, this would give (Eq. [T0| 

a;^(d = 2) = a;(d = 2) = i (32) 

in agreement with the measure of the critical typical correlation via Monte-Carlo [l5l |. 

As a final remark, let us stress that even in dimension d — 1, the critical behavior of the averaged correlation 

Ccriti{r) OC r"^^ (33) 

r— f +00 

is given by a non-trivial power-law (involving '^-[r^'jj^ = Pid = 1) = (3 — \/5)/2) that cannot be obtained simply 
from the Brownian scaling. This is because it represents some 'persistence exponent' [2^ for the strong disorder 
renormalization, that can only be computed from the exact solution of the renormalization flow . As a consequence, 
we do not expect any simple exponent for the critical averaged correlation in dimension d > 1 either. 

IV. STATISTICS OF THE LOCAL SUSCEPTIBILITY WHEN aj(d) > 
A. Perturbative expansion deep in the disordered phase of the quantum model 

In the presence of a uniform infinitesimal exterior field id — 0, the local magnetizations rrii ~ xioc{i)H at linear 
order in H defines the local susceptibilities xioc{i)- At first order in the couplings J,;.j, the local susceptibility Xioc(O) 
is related to the local susceptibilities xioc{j) of its K = {2d) neighbors on the hypercubic lattice via 

1 T 

Xioc{0) = j^+Y.f^Xioc{j) (34) 
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The difference with Eq |3] is the first term I/Hq that represents the local susceptibility of the isolated spin (when 
Ji,j =0)- If we iterate Eq[Ml we obtain 

xioc =J2 E ^(^) 

r PathP-.O^r 

Win ={ U ^] (35) 

The difference with Eq|3]is that the path starting at point now ends at any point r in the volume L^. The points r 
very close to give finite contributions, whereas the points f far away will have a contribution displaying the scaling 
(see Eqs Hand III) 

J2 W{P) ~ 6"^+"^"^^ (36) 

PathP-.O^r 

where the random variable u — up is of order 0(1). 



B. Power-law distribution of the local susceptibility 



A saddle point analysis of Eq. [36] yields that the leading contribution will come from points r at a distance 

r saddle ^ iS.typU)T~^ = S.avU^ (37) 

i.e. of the order of the averaged correlation length £^av of Eq. [26l This saddle point gives a contribution of order 

In X;oc - + ^l^ur^ = (38) 
It is thus useful to introduce the positive random variable 

_ In Xioc 1 

V = —— — oc u^-^ — u' (39) 

From the tail of the distribution of the variable u of Eq. [5T1 one obtains that the distribution Q{v) of the variable 
V cx u'' presents the exponential tail (Eq. [5T] ) 

Q{v) cx e""" (40) 

where a is some constant. The change of variables from v to xioc of Eq. [39] yields that the distribution of the local 
susceptibility presents the power-law tail 

Pixioc) 'x (41) 

where the exponent /i varies continuously with the averaged correlation length ^av 

M = 7^ (42) 

As the critical point is approached ^av +oo, this exponent /i goes to zero, and the distribution P{xioc) becomes 
broader and broader. In particular, in the region near the critical point where < /i < 1, the averaged local 
susceptibility diverges xzoc = +oo (Griffiths phase), whereas farther away where u > 1, the averaged value xi 
remains finite. These conclusions again generalize directly what is known in d = 1 
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C. Scaling of the local magnetization as a function of the magnetic field 

In the Griffiths phase < /i < 1 where the averaged locai susceptibiiity diverges x/oc — +oo, this means that the 
averaged magnetization rn wiil grow more rapidiy than linearly in the magnetic field H. To see what happens, it is 
convenient to consider the probability to have a magnetization mg = XiocH of order 1 from Eq. 21] 

Prob{mo = xiocH - 1) = ^Prob{xioc > ^) (43) 
The averaged value WTo will be thus governed by these rare events and scales as 

rn^ - cx Hf" ^ (44) 



In summary, the analysis of this section generalizes the exact results in c? = 1 Q to all finite dimensions d where the 
droplet exponent remains finite io{d) > 0, in order to ensure the vanishing of the exponent fi of Eq. H^lnear criticality. 



V. DYNAMICAL SCALING WHEN > 



An essential property of 'conventional' quantum systems is the dynamical exponent z that describes how the gap G 
(that defines the appropriate characteristic inverse time) of a sample of volume L'^ vanishes with L as the power-law 

G(L'^) - L^' (45) 

For instance the pure quantum Ising model of Eq. [T]is characterized by Zp^re = 1 in any dimension d, as a consequence 
of the equivalence with a classical Ising model in dimension (d + 1), where the time plays the role of an extra spatial 
dimension [1,]. For the disordered model of Eq. [5] in dimension d = 1, the exact solution [^f corresponds however to 
the following activated scaling at criticality 

lnG(i'^) - -L'^^'*) (46) 

with •ip{d = 1) = 1/2, instead of the power-law of Eq. 1451 meaning that the dynamical exponent z is formally infinite 
at criticality. In the following, we thus discuss the scaling properties of the gap in the disordered phase in dimension 
d > 1 where w(d) > 0. 



A. Gap of a sample of volume L'' in the disordered phase 



To relate the study of the gap with the study of the local susceptibility discussed in the previous section, it is 
convenient to follow Ref [s^ and to write the local susceptibility in terms of the many-body eigenstates |n > of 
energies i?„, where n = labels the ground state (Eq. 67 of Ref |30| ) 

«-« = 2i:li|Kti£ ,47) 

n>0 " ° 

One expects that xioc will behave like the inverse local gap Gioc, i-e- more precisely, that the local gap Gioc in a 
region of volume (Eq. [ST]) 

V ^ rtddie ~ ?i (48) 

behaves as 

Gioc - — (49) 

Xloc 

The power-law tail of the distribution P{xioc) of Eq. 21] transforms into the following singularity for the distribution 
P{Gioc) at the origin 

P(Gloc)^^^^A^^G';-^' (50) 
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where A is some amplitude. Note in particular that in the Griffiths phase < /i < 1, this correspond to a diverging 
singularity at the origin. 

As discussed in detail in [bU. l32j|. the gap G{L'^) of the whole system of volume L'^ in the disordered phase is given 
by the minimum value G]^™ of the N — L'^/^av independent local gaps Gioc that are present in the system. Taking 
into account the power-law singularity of Eq. [50l one obtains that the gap G{L'^) = Gmin scales as 



1 



dGiocPiGioc)- A{GrorT (51) 



i.e. by inversion 







GiL'') = GIT ^A-(tj] (52) 



From the dependence in L (Eg I45p . one reads that the dynamical exponent z is directly related to the exponent fi 
discussed in the previous section (Eq. HT|) 

z = - (53) 

This relation involving the factor d [H, [l^ generalizes the well-known result z — l/fi in d = 1 [2j]. In particular the 
Griffiths phase < /i < 1 corresponds to the domain z > d. The vanishing of the exponent ^ near criticality (Eq 
yields the following divergence of the dynamical exponent near the transition 

- = - = -C (54) 
fi a 

B. Finite-size scaling for the typical gap 

From the divergence of the dynamical exponent in Eq. I54[ we expect that the critical point cannot be conventional 
but should obey the activated scaling of Eq. |46]with some exponent ip. In addition, all finite-size-scaling properties 
are expected to be governed by the averaged correlation length £^av , as in Eq. [28l for the typical correlation. It is thus 
natural to assume the following finite-size scaling form for the typical gap, in agreement with the exact finite-size 
scaling form known in d = 1 [3l| 

L 



\nG{Ld) ^ ^L^<i>l \ (55) 

where the scaling function reduces to some positive constant $(0) > at criticality. The matching with the scaling 
in the disordered phase of Eq. [52] (using the critical behavior of the exponent of Eq. 



(I) (I 



(56) 



imposes the asymptotic behavior 



$(a;) oc x-'^ ln(a;'') (57) 



leading to the conclusion (see also lSTj) 



1p = U! = iUc (58) 

that generalizes the known result ip = oj = 1/2 in d = 1 2]. In dimension d = 2, this would give (Eq. [TU)) 

ij{d = 2) = uj{d = 2) = \ (59) 

The numerical estimations of ip tend to be larger, around ip ~ 0.42 [1, [l^ or ^/i ~ 0.48 [ll|, whereas the measure 
of U!c of typical critical correlation, which is supposed to be the same fEg [5T|) is closer to 1/3 [lH|. Further work is 
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needed to better understand the origin of these differences in d = 2, as well as in d = 3,4 where the measures of 
the gap exponent tp [H, [13] remain around 0.5, whereas the droplet exponent oj{d) = ujdp{D = d — 1) decays with 
d (see Eq. lASp . A possible explanation for this discrepancy could be that all these Strong Disorder RG numerical 
results are based on the approximation of the 'maximum rule' introduced in the very first paper concerning d > 1 
that allows to study much bigger sizes. This 'maximum rule' approximation has been recently questionned [33| for 
another quantum model, namely the superfluid transition for random bosons (see the discussions in section III B and 
Appendix A of [33|). For our present model, we believe that the use of the approximated 'maximum rules' instead of 
the full 'sum rules' could be a problem for the following reasons : 

(i) the full 'sum rule' corresponds for the Directed Polymer model to the computation of the partition function as 
a sum over all paths, so that what shows up is the fluctuation exponent Ld^p < 1/2 of the free-energy, that we have 
discussed in detail in the present paper. 

(ii) the approximated 'maximum rule' corresponds for the Directed Polymer model to the computation of the 
contribution of a single dominant path, i.e. of the energy of this path. However at all non-zero temperature T > 0, 
the energy and the entropy of the Directed Polymer are known to display fluctuations of order L^/^ in all dimensions d, 
coming for independent short-scale contributions. There exists a subtle 'cancellation' between the fluctuating parts of 
the energy and the entropy at leading order, so that the free-energy displays smaller fluctuations of order L'^op ^ 2^1/2 
(49I that leads to the phenomenon of 'chaos' in temperature [i^, as in spin-glasses [soj . 

To confirm this interpretation, one needs to study numerically the full 'sum rule', but unfortunately this seems to 
be possible only for too small system sizes to get reliable estimations of the exponents. The remaining possibility 
could be to develop other types of approximation, like the cut-off approximation introduced recently in [33|. 



VI. CASE OF THE CAYLEY TREE 



As already mentioned in section lll El the physics of the random field Ising model on the Cay ley tree in its disordered 
phase has been related to the physics of the Directed Polymer on the Cayley tree in Refs [23J25| via the quantum cavity 
method. In this section, we mention the consequences of this correspondence for various observables to emphasize the 
similarities and differences with the case of finite dimension d. To allow for the possibility of the two phases (delocalized 
or localized) for the effective Directed Polymer, we have chosen here not to consider the 'Box'-distribution of Eq. [11] 
(which is discussed in detail in Refs ^23n25j ). but to focus instead on the following log-normal distribution of the 
random fields 

1 (Inh-hThj^ 

^LNih) = e (60) 

AlV 27rCT^ 

of fixed parameters {hih,a), whereas the ferromagnetic couplings Jij are not random but take a single value J that 
will be the control parameter of the quantum transition. 



A. Magnetization of a center site 

On a Cayley tree of coordination number (K + 1) with L generations, the local magnetization satisfies the following 
recurrence at lowest order in perturbation theory (analog of Eq. [3]) 

J 

rriL («) = 7- X! (61) 

' i=i 

where m^fLi are K independent realizations of the magnetization after [L — 1) generations. This corresponds to the 
recurrence of the Directed Polymer on a Cayley tree (Eq. IA15[) with temperature T = 1 and with effective random 
energies (instead of Eqs H] and 

e, = In /i, - In J (62) 

Eq. ISni corresponds to the Gaussian distribution of Eq. IA14I of averaged value 

eo = In /i — In J (63) 

The translation of the exact Derrida-Spohn solution [4&| recalled in section IA4I of the Appendix, yields that there 
exists a critical width Oc for the disorder distribution of Eq. [60] where the critical temperature of Eq. IA19I is unity, 
i.e. 



CTc = y/2\a.K 



(64) 



12 



1. When the Directed Polymer is in its delocalized phase a < 

The Directed Polymer is in its delocalized phase atT = lifT=l> Tc, i.e. below the critical disorder a < ac 
Eq. R20I becomes at T = 1 



In if + — + In J - hiTi 



i + M (65) 

K 'H-ext / -i^— >+oo I z J 

where the distribution P'^'^^°'^{u) of the random variable u presents the exponential tail (Eq. IA2ip 

P'^'='°=(u) cx e~(^)" (66) 
Equivalently, if one introduces the typical value (Eq. [65] ) 

typ _ Inmn -(\ah-\aK-^~\nj\L 

ttLq = e ' CX rriext e V 2 y (^gyj 

one obtains that the distribution of the magnetization too presents the following power-law tail 



P'^«'-(too) ^ CX, ^ with Trfeio, = ^ > 1 (68) 



From Tdeioc > 1, one obtains that the ratio toq/toq''^ remains finite, so that the averaged value has the same behavior 
as the typical value of Eq. IHT] . 

If one assumes that this Directed Polymer approximation remains valid in the whole disordered phase of the quantum 
model, the behavior of the typical magnetization of Eq. [57] yields that the critical coupling Jj, would be 

2 

In Jc = h^ - In if - y (69) 

and that the radial correlation length describing the exponential decay in the radial direction of the tree would be 
given for J < Jc by 

= In Jc - In J cx {Jc - J)"'- with i/^ = 1 (70) 

2. When the Directed Polymer is in its localized phase a > ac 

The Directed Polymer is in its Localized Phase at T = 1 if T = 1 < Tc, i.e. above the critical disorder cr > Uc. Eq. 
E22I becomes at T = 1 

In ( ) ~ InZDp(L) =i \cjC7c + In J -]nh\ L- —In L + u (71) 

\mext J 2ac 

where the distribution P^°^{u) of the random variable u presents the exponential tail (Eg IA23|) 

F'°"(u) cx we""^" (72) 

Equivalently, if one introduces the typical value (Eq. [7T|) 

TO*^P = e'^ oc TOe^t L-*^ e-('^^""=-^"'^)^ (73) 

L— )- + oo 

one obtains that the distribution of the magnetization toq presents the following power-law tail (with a logarithmic 
correction) 

V'°'{m,) « ijH^Tz 1, ^ith = ^ < 1 (74) 
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From Tioc < 1, one obtains that the ratio toq jniQ^ diverges. 

If one assumes that this Directed Polymer approximation remains valid in the whole disordered phase of the quantum 
model, the typical magnetization of Eq. [73] yields the critical coupling Jc would be 

In Jc = In /i — CTCTc (75) 

and that the radial correlation length for J < Jc would be 

— = In Jc — In J cx ( Jc — JY^ with Vr = 1 (76) 

3. Possible finite-size scaling 

Critical points of disordered models on the Cayley tree usually involve two correlation leiigth exponents, as shown 
by the examples of the localization/delocalization transition of the Directed Polymer [sj, [1^, the random wetting 
transition [sO] and of Anderson transition [37} . From the analogy with the Anderson transition, where the localized 
phase can be studied by a traveling-wave approach, and where the transition corresponds to a traveling-non traveling 
phase transition, one might expect the following finite-size scaling form for the typical center-site magnetization 



In ( ^ ) ~ -LP^ (l^/^^^Jc - J)) (77) 

yrriext J L^+oo \ / 

with some finite-size scaling exponent i'fs- This type of finite-size scaling for a traveling/ non-traveling phase transition 
has been found in Ref 38] for a soluble case with explicit critical exponents, and for the Anderson transition with 



numerical estimates of the exponents 37[. The matching with the disordered phase involving the radial correlation 



length ( Jc — J) discussed previously 



leads to the relation 



In(^) c -L^^LiJc~Jr (78) 



p+^ = l (79) 



In the ordered phase J > Jc, one might expect in analogy with [s^, the essential singularity 



In ^ -(J-Jc)-- (80) 

The matching with Eq. [77] gives the relation 

K = PI'FS (81) 

Within our present approach, we can use the associated Directed Polymer model only in the quantum disordered 
phase, but not at criticality. As a consequence, we cannot justify the behavior at criticality or the essential singularity 
of Eq. 1801 Using some approximations within the quantum cavity method (which may however break down close to 
criticality), Refs (23l - l25j have found the essential singularity of Eq. [5n]with k = 1. Taking into account Vr = 1 (Eqs 
[75|) . this would correspond to the simple values vps — 2 and p = 1/2. 

As a final remark, we note that for the Anderson transition, the exponent k = 1/2 of the essential singularity has 
been exactly computed (39j 

B. Statistics of the local susceptibility in the disordered phase 

On the Cayley tree, the recurrence at lowest order for the local susceptibilities read (Eq. [34| 

^ioiW = )!- + £Exir^W (82) 
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where Xioc~^Hj) ^re K independent realizations of the local susceptibilities after (L — 1) generations. Without 
the inhomogeneous term l/h^, one would recover the recurrence of Eq. [HI] for the magnetizations that decay 
exponentially towards zero in the disordered phase. As a consequence, one expects that the recurrence of Eq. [5^ 
converges towards a finite random variable xioc = Xio^'^°°'' that should be stable upon iteration. Upon iteration, Eq. 
[83l yields that the local susceptibility can be seen of a sum over partition functions of the Directed Polymer on the 
tree of arbitrary length r (analog of Eq [35]) 

Xioc-Y.^Dpir) (83) 

r 

Let us discuss the behavior of the averaged value xum, the behavior of the typical value xloc^ ^^'^ then the tail of 
the probability distribution P{xioc)- 

1. Behavior of the averaged value 

Averaging Eq. [82] over the disorder leads to the following closed expression for the averaged value xioc, that we 
may explicitly compute using the disorder distribution of Eq. 1601 and the value of Jc of Eq. 

Jdh^^^^ _ e^-i^ _ 1 

i.e. it diverges with the 'mean- field' exponent 7 = 1. 

2. Typical value for a < 

For a < <Jc where Eq. 1551 holds 

ZDP{r) ~ e-*+" (85) 
we expect that the typical local correlation scales as the radial correlation length 

~ / drZopir) ^ I dre-^ cx oc ( Je - J)-"'-" (86) 



Taking into account vtyp = 1, one obtains that the typical value diverges with the same mean- field exponent jtyp = 1 
as the averaged value of Eq. [M] 



For a > (Tc where Eq. [TT] holds 



one obtains the following contribution 



Typical value for a > ac 



ZDp{r)~r-—^ e--+" (87) 



X;oc J drZopir) ^ J drr e oc 

that does not diverge with since a > CTc- This means that the typical value Xioc remains finite as the transition 
is approached +oo, as a consequence of the value of the coefficient > 1 of the logarithmic correction of the 
Directed Polymer model (Eq. IA22I) ! 
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4- Exponent /j, of the power-law distribution 



From the recurrence of Eq. [SU one expect that the stable distribution P{xioc) will display a power-law tail 



PiXloc) 



1+p 



(89) 



where A is some amplitude. Note that for the special case K = 1 where the tree degenerates into a single line, the 
recurrence of Eq. |82] defines a so-called Kesten variable that has been much studied in various physical contexts 
[40| - |43 | with the conclusion that /i is determined by the condition (written for our present case) 



Case K = 1 



(90) 



For the tree with ii' > 1, it is straightforward to adapt the argument as follows. The distribution P{xioc) has to be 
stable upon the iteration of Eq. 1821 when the random field hi is distributed with i^LNih) 



K 



(91) 



The stability of the power-law tail of Eq |89] in the region xio 



-oo can be analyzed as follows : a very large xio 



is obtained only if one of the K values Xiic large; assuming it is = 1, one has then xio 

function, so that Eq IHI] becomes in the tail region 

K ldhnMh)ldxj2^^^s{x..-'j^xt 



J U) 

h^loc 



in the delta 



A{xT^^ 



Xi 



^loc ' 



J dhTTLNih)h- 



(xiocy 

{XlocV+^ 



yielding that the tail exponent n is stable only if it satisfies 



that generalizes Eq. |90]to arbitrary K. 

With the disorder distribution of Eq. [BDl the stability condition of Eq. 



1 reads 



1 — ^g^^^-A'(lii'»-ln-f) 

Taking the appropriate solution of this second-order equation in /i, one finally obtains 



1 



(In/i-ln J) + J(ln/i-ln J)2 -2cr2 InK 



(92) 



(93) 



(94) 



(95) 



that generalizes the known result /i = ^ (In — In J) for the one-dimensional case K — 1 Kesten variables. To see 
more clearly how the exponent /i varies within the quantum disordered phase as a function of the radial correlation 
length S.r, let us now distinguish the two cases : 
(i) For a < ac, using Eqs [64l and iTOl one obtains 



(96) 



The exponent fideioc decreases as grows and the limit value near criticality +oo 



16 



coincides with the exponent Tdeioc > 1 of Eq. 1681 
(ii) For (T > (Tc using EgslMlandlTSl 



_ 1 




The exponent /ijoc decreases as grows and the limit value near criticality 



-oo 



^ 1 



1l 

a 



(98) 



(99) 



coincides with the exponent tioc < 1 of Eq. [731 

In summary, for the log- normal distribution of the random fields (Eq. [60]) . the statistics of the local susceptibility 
for (T > (Tc is very similar to the results found in [23l - [25j for the case of the Box Distribution of Eq. [TT] : the exponent 
Hioc converges towards a finite value — < 1 as the critical point is approached — > +oo, instead of converging 
towards zero when Lj{d) > in finite dimension d (Eq. I42|) . 



C. Dynamical scaling in the disordered phase 

As in section |Vl the power-law distribution of the local susceptibility (Eq. [89)) leads to a power-law singularity for 
the distribution P{Gioc) of the local gap Gioc (see the discussion leading to Eq. [33]) 

P{Gioc) ^^^^A{GZ^^)-^GI-' (100) 

In a Cay ley tree with L generations of volume , the number N of independent local gaps Gioc will also grow 
exponentially in L. As a consequence, the gap G{K^) of the whole system in the quantum disordered phase that is 
given by the minimum value GJ^"" of the N independent local gaps Gioc is expected to scale as 

InG(X^) = lnG[2™ cx In f cx --IniV cx ~L (101) 

instead of Eq. [521 in finite dimension d. So the gap G{K^) decays only exponentially with respect to the radial size 
L, and this corresponds formally to an activated exponent -0 = 1 and to an infinite dynamical exponent z = -\-oo 
already within the quantum disordered phase, in agreement with the limit d — > -|-oo of Eq. [S31at fixed /i. 



VII. CONCLUSION 



For the random transverse field Ising model in finite dimensions d > 1, we have analyzed the lowest-order con- 
tributions in perturbation theory in ( Jij/hi) to obtain some information on the statistics of various observables in 
the disordered phase. We have established a link with the statistics of the partition function of the Directed Poly- 
mer with D = (d — 1) transverse directions. We have analyzed the consequences in terms of the droplet exponent 
uj — ijj£,p{D — d — \) oi the associated Directed Polymer with the following conclusions : 

(i) Whenever the Directed Polymer is in its localized phase with a positive droplet exponent a; > 0, the quantum 
model is governed by an Infinite-Disorder fixed point. In particular, there are two distinct correlation length exponents 
related by vtyp = (1 — ui)i'av The distribution of the local susceptibility xioc presents the power-law tail P{xioc) ^ 
1 /Xi^c^ where vanishes as , so that the averaged local susceptibility diverges in a finite neighborhood < /i < 1 
before criticality (Griffiths phase). The dynamical exponent z diverges near criticality as z = d/ ii ^ 

(ii) In dimensions d < 3, the associated Directed Polymer is always in its localized phase (the delocalized phase 
does not exist), so that any infinitesimal disorder flows towards this Infinite-Disorder fixed point with Lu{d) > (for 
instance Lu{d = 2) = 1/3 and uj{d = 3) 0.24) 

(iii) In finite dimensions d > 3, the associated Directed Polymer can be in two phases depending on the disorder 
strength : 

(iiia) For strong enough disorder, the fiow is towards an Infinite-Disorder fixed point provided the droplet exponent 
of the localized phase of the associated Directed Polymer remains positive lu = ujdp{D = d — 1) > (this seems the 
case for at least d < 8 from the numerical results quoted in Eq. IA5I and it could be for all d < +oo if the upper 
critical dimension is infinity for the Directed Polymer model). 
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(iiib) For small enough disorder, the associated Directed Polymer will be in its delocalized phase with utdeioc = 0, 
so that the quantum model could then flow towards a more conventional Finite-Disorder fixed point. 



Note added in proof : 

In this paper, we have presented the simplest scenario where the random variable u of Eq. 7 or 20 remains of order 
0(1) as in the exact solution in dimension d=l and as found in dimension d=2 via a non-linear transfert approach 
[C. Monthus ans T. Garel, J. Stat. Mech. (2012) P01008]. However, a more complicated scenario where the random 
variable u of Eq 7 or 20 contains a diverging amplitude as the critical point is approached, is also possible, as was 
found recentlty on a hierarchical fractal lattice [C. Monthus ans T. Garel, arXiv:1201.6136 . We refer to this preprint 
for the description of the consequences on various scalings and on relations between critical exponents. 



Appendix A: Reminder on the Directed Polymer in a random medium 

The model of the Directed Polymer in a (1-1- D) random medium (where D is the dimension of transverse spatial 
directions, and where the 1 refers to the 'directed' direction) is defined by the following classical partition function at 
inverse temperature /? = 1 /T 

^r(/3) = E^^p(-/3 E (Ai) 

RW \ l<a<L J 

The sum is over D— dimensional random walks f{a), where the effective 'time' a represents the directed direction. 
The independent random energies e(a, r) can be put either on the sites or on the bounds. This model has attracted a 
lot of attention because it is directly related to non-equilibrium properties of growth models in the KPZ universality 
class [45|. Within the field of disordered systems, it is also very interesting on its own because it represents a 
'baby-spin-glass' model [isl - ligj . 

The possible phases of the Directed Polymer are : 

(i) a localized phase at low temperature T < T^. where the order parameter is an 'overlap' [1^, H^. In finite 
dimensions, a scaling droplet theory was proposed [49|, in direct correspondence with the droplet theory of spin- 
glasses [131, whereas in the mean-field version of the model on the Cayley, the localized phase is very similar to the 
frozen phase occurring in the Random Energy Model (46j . 

(ii) a delocalized phase at high temperature T > 

In the following, we recall the properties of these phases that are needed in the text, as well as the phase diagram 
as a function of the dimension D 



1. Properties of the localized phase 

a. Droplet exponent u}dp{D) 

The droplet theory for Directed Polymers 49] is similar to the droplet theory of spin-glasses [s^l- It is a scaling 
theory that can be summarized as follows. At very low temperature T — > 0, all observables are governed by the 
statistics of low energy excitations above the ground state. An excitation of large length / costs a random energy 

Ai;(0 - (A2) 

The droplet exponent w^ip is exactly known in dimension D = 1 [ssl [sgI Is^ . l65j 

ujDp{D = 1) = ^ (A3) 
and for the mean-field version on the Cayley tree [i^ 



(A4) 
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In finite dimensions D = 2,3,4,5, the exponent uj]jp{D) has been numerically measured, in particular [57H6i 



ujdp{D = 


2) 


~ 0.244 


(jJDp{D = 


3) 


~ 0.186 




4) 


~ 0.153 




5) 


~ 0.130 


ujdp{D = 


6) 


~ 0.114 


^Dp{D = 


7) 


~ 0.100 



(A5) 



The authors of Ref. [61| have moreover argued that the droplet exponent remains positive ujdp{D) > in all finite 
dimension D. 



b. Distribution of the partition function : tail exponent rj 

One expects that the whole low temperature phase < T < Tc is governed by the zero-temperature fixed point 
characterized by the droplet exponent ujdp- In particular, the logarithm of the partition function scales for T < Tc as 

IuZl ~ -pfiocL + L'^^^u (A6) 

The first term is extensive and non-random {fioc represents the free-energy per step in this localized phase). The 
second term involving the droplet exponent lodp is random : u is a random variable distributed with some law 
that presents the following asymptotic behavior 

Inp(u) cx -u^ (A7) 



where the tail exponent ry(-D) is directly related to the droplet exponent ujp,p{D) via (see [62, [ZO, IZll and references 
therein) 

VDP = --^ (A8) 

1 - UJDP 



In dimension _D = 1, the distribution p{u) is exactly known to be the Tracy- Widom distribution [6J-|69| with the 
tail exponent 

r,DP{D = l) = l (A9) 
On the Cayley tree, the distribution p{u) is also exactly known |46,] with the tail exponent 

v'£F = 1 (Aio) 



2. Properties of the delocalized phase 



In the delocalized phase T > Tc , the logarithm of the partition function scales as (compare with Eq. IA6p 

IuZl ^ -PfdelocL + U (All) 

In the extensive non-random term, the free-energy per step fdeioc is given by the annealed value fann that can be 
obtained from the averaged partition function Zl 



( HZl) 



P fdeioc = -13 fann = Hm (A12) 



The second term in Eq. lAlll is a random variable u of order L°. 
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3. Possible phases as a function of the number D of transverse directions 

The phase diagram as a function of space dimension D is the following (45j : 

(i) In dimension D < 2, there is no delocalized phase i.e. the critical temperature is infinite 

Tc{D < 2) +00 (A13) 

and any initial disorder drives the polymer into the localized phase. 

(ii) In dimension D > 2, there can be a phase transition between the low temperature localized phase and a 
delocalized phase at high temperature [sil, H^l ■ This phase transition has been studied exactly on a Cayley tree . 
In finite dimensions, bounds on the critical temperature Tc have been derived [53 - l&3 | : Tq < Tc < T2. The upper 
bound T2 corresponds to the temperature above which the ratio remains finite as L — t- 00. The lower 
bound Tq corresponds to the temperature below which the annealed entropy becomes negative. For instance, when 
the random energies in Eq. lAll are Gaussian variables 

the temperature T2 can be exactly computed and is finite for D > 5 [53J , so that one is sure that the delocalized phase 
exists at least in the region T > 



4. Reminder on the exact solution on the Cayley tree [4] 

The Directed Polymer model on the Cayley tree of coordination number (K + 1) can be solved for any disorder 
distribution via the Derrida-Spohn traveling- wave approach (46| for the recurrence 

K 

Zi(z) = e-^^-^Zi_i(j) (A15) 

where z'j^}_^ are K independent realizations of the partition function after (L — 1) generations. Here we recall the 
solution for the Gaussian distribution of Eq. IA14l that leads to very simple explicit results : the annealed partition 
function reads 

Z^= (^K^y = (^Ke^i^"")^ (A16) 
so that the annealed free-energy per step at temperature T = 1//3 reads 

UUT) ^ ^Hm^ (-^) - e„ - TlnX - ^ (A17) 
with the corresponding annealed entropy 

2 

SanniT) = -drfanniT) = \n K - ^ (A18) 

The transition temperature Tc corresponds to the temperature where the annealed entropy vanishes Sa„„(Tc) = so 

Tc - , " (A19) 

and the two phases have the following properties 

In the Delocalized Phase T > Tc, the free-energy coincides with the annealed free-energy of Eq. IA171 and the 
logarithm of the partition function has for statistics 



In^L ^ 

L— >+oo 



L + u (A20) 
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where the random variable u of order 0(1) has a distribution P"^'°'^(w) that presents the exponential tail [4] 



P'"'""{u) oc e V^. 



(A21) 



In the Localized Phase T < Tc, the free-energy per step remains frozen at the value fann{Tc), and the logarithm of 
the partition function has for statistics 



InZf 



L— >+oo 



f T 



ST 

L ■ InL + u 

2T 



(A22) 



where the first correction to the extensive term is a non-random logarithmic term |34l . |46j , and where the variable u 
of order 0(1) has distribution P^°'^{u) that presents the exponential tail ^1] 



u— >-+oo 



(A23) 



Appendix B: Analogy with the localized phase of Anderson localization 

The Anderson tight-binding model for a single quantum particle on a hypercubic lattice in dimension d is defined 
in terms of the Hamiltonian 

H = Y,e.,\z><z\+Y,V,,jV><j\ (Bl) 

where ti is the random on-site energy on site i and where V^ j- is the hopping between the sites i and j (usually taken 
to be unity V = 1 between nearest neighbors). 

In dimension d = 1, the transfer matrix formulation of the Schrodinger equation yields a log- normal distribution 
for the Landauer transmission Tl [zl, 

lnTl''=^> oc -A + ii/2„ (B2) 

The leading non-random term is extensive in L and involves the localization length ^/oc- The subleading random term 
is of order L^/^, and the random variable u of order 0(1) is Gaussian distributed as a consequence of the Central 
Limit theorem. Although it has been very often assumed and written that this log-normal distribution per sists in the 
localized phase in dimension d = 2,3, theoretical arguments [t^, [lE] and recent numerical calculations |76l. are in 
favor of the following scaling form for the logarithm of the transmission 

\yiT^^'^ cx (B3) 

where the exponent uj{d) depends on the dimension d and coincides with the droplet exponent characterizing the strong 
disorder phase of the directed polymer in a random medium of dimension \ + D with D — d — 1 (see the reminder 
on Directed Polymers in Appendix |^ . The probability distribution of the rescaled variable u is not Gaussian but is 
determined by the directed polymer universality class (see [76| where its distribution in d = 2 is shown to coincide 
with the exactly known Tracy- Widom distribution for the directed polymer in 1 -I- 1). 

The arguments in favor of the same universality class can be decomposed in two steps [T^^ZHI : 

(i) in the localized phase of Anderson localization in dimension d, the transmission decays exponentially with the 
length, and thus directed paths completely dominate asymptotically over non-directed paths. 

(ii) these directed paths of the Anderson model have weights that are random both in magnitude and sign, but it 
turns out that the directed polymer model which is usually defined with random positive weights (Boltzmann weights) 
keeps the same exponents in the presence of complex weights (see section 6.3 of the review |45|). 

Note that the Directed Polymer model fixes the fluctuation exponent w, the rescaled distribution of the random 
variable u, but does not give information on the divergence of the localization length S^ioc near the Anderson transition. 
We expect that the same conclusions hold for the random transverse fields Ising model (see section HTCj) . 
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